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Abstract 


Heavy metals are known to be the most dangerous environmental pollutants because they do not 
decomposed by physical processes and therefore remain for a long time and will be effective in 
biochemical cycles and ultimately in the human food chain, disrupting biological reactions and 
damaging organs. Therefore, in order to control the quality of the environment, it is important to 
study the heavy elements in the soil. So far, many researchs have been done to identify 
contaminated areas and implement sustainable land use policies regarding the risk of 
contamination. In this study, the spatial distribution of nickel in surface soils (0-10 cm) in parts 
of Qom province with urban-industrial, agricultural and uncultivated land uses with an area of 
883 km? was investigated. Geostatistical tools (ordinary kriging) were used. Sampling was 
performed at 209 points on a grid with distances of about 1.5 x 1.5 km in agricultural and urban 
land uses and about 2x2 km in uncultivated lands. The concentration of total nickel and 
absorbable soil was measured by atomic absorption spectrometer. Statistical and geostatistical 
measures were performed by SPSS, Variowin and WINGSLIB software and contamination map 
was drawn with Surfer16 software. The results indicated that the concentration of nickel in various 
forms in the region does not show toxicity. But there were significant differences in land uses. 
Existing industries, parent materials, and nickel-rich basic and ultra-basic rocks, the use of 
phosphate fertilizers and sewage sludge may have increased nickel, and the direction of westerly 
winds has led to the transfer of nickel in the region. 
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جغرافیا و مخاطرات محیطی. سال یازدهم. شمارة چهل و یکم بهار ۶۰۱ صص ۱۷-۳۷ 
مقاله پژوهشی 


توزیع منطقه‌ای نیکل در خاک‌های سطحی بخش‌هایی از استان قم 


زهرا موحدی راد - مربی گروه کشاورزی دانشگاه OLS og ely‏ ایران 
(IB)‏ سمیه صدر' - استادیار گروه کشاورزی؛ دانشگاه ply‏ نو LS‏ ایران 


تاریخ دریافت: ۱۶۰۰/۲/۸ تاریخ بازنگری: ۱۶۰۰/۳/۱۱ تاریخ تصویب: ۱۶۰۰/۶/۱۸ 


ce 


تا کنون تحقیقات بسیاری در زمینه شناسایی مناطق آلوده و اجرای سیاست‌های استفاده پایدار از اراضی با 
توجه به خطر آلودگی آن‌ها انجام گرفته است. در این پژوهش توزیع مکانی نیکل در خاک‌های سطحی 
Gas)‏ ۰-۱۰ سانتیمتر) بخشی از استان قم با کاربری‌های شهری-صنعتیء کشاورزی و بایر به وسعت MÝ‏ 
کیلومترمربع بررسی شد. در این مطالعه از ابزار زمین‌آمار (کریجینگ معمولی) استفاده گردید. نمونه‌برداری 
در ۲۰۹ نقطه بر روی شبکه‌ای با فواصل حدود Vox MO‏ کیلومتر در اراضی کشاورزی و شهری-صنعتی و 
حدود ۲×۲ کیلومتر در اراضی بایر انجام شد. غلظت نیکل کل و قابل جذب SLE‏ توسط دستگاه جذب 
اتمی ان دازه‌گیری شد. مطالعات آماری و زمینآماری توسط نرم‌افزارهای Variowin SPSS‏ و 
WINGSLIB‏ انجام و نقشه آلودگی با نرم‌افزار 5۱116716 رسم گردید. نتایج نشان داد که غلظت عنصر 
در هیچ‌یک از شکل‌ها در منطقه سمیت را OUS‏ نداد؛ اما در کاربری‌های مختلف دارای تفاوت معنی‌دار 
بود. احتمالاً صنایع موجود. مواد مادری و سنگ‌های بازیک و فوق بازیک که غنی از نیکل هستند. استفاده 
از کودهای فسفاته و لجن فاضلاب باعث افزایش JSS‏ و جهت وزش بادهای غربی باعث انتقال نیکل در 
منطقه گردید است. 
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۱- مقدمه 

سلامت محیط‌زیست موجودات زنده به خصوصیات شیمیایی خاک به‌حصوص polite‏ کافی عناصر سنگین 
ارتباط نزدیک دارد. امروزه در بسیاری از کشورهاء اهداف جدید قانونگذاری» ارزیابی و کنترل فلزات سنگین در 
خاک است؛ چراکه در دهه‌های اخیر مسئله آلوده شدن خاک‌ها به عناصر سنگین در سطح وسیع» به یک امر مهم 
زیست‌محیطی تبدیل شده است (اسلام ! و همکاران, ۲۰۱۵). 

فلزات سنگین به‌عنوان حطرناک‌ترین آلاینده‌های زیست‌محیطی شناخته می‌شوند؛ چراکه به علت عدم تجزیه 
توسط فرایندهای فیزیکی, به مدت طولانی در محیط باقی‌مانده و در چرخه‌های بیوشیمیایی و نهایتاً زنجیره غذای 
انسان موثر خواهند بود و سبب اختلال در واکنش‌های بیولوژیکی. آسیب به اندام‌ها و حتی مرگ خواهند شد 
(کاساسی " و همکاران» ۲۰۰۸)؛ بنابراین جهت کنترل کیفیت محیط‌زیست بررسی عناصر سنگین در SLE‏ دارای 
اهمیت است. اگرچه برخی از عناصر سنگین به‌عنوان عناصر غذایی ضروری هستند. اما اغلب آن‌ها به علت غلظت 
UL‏ در خاک و یا پتانسیل قابلیت دسترسی زیستیء برای بیوسفر سمیت داشته و خطرناک می‌باشند (کایلود "و 
همکاران, ۲۰۰۹؛ ساری* و همکاران: ۲۰۱۸). ازجمله این عناصر می‌توان به نیکل اشاره کرد. این عنصر به‌واسطه 
ناراحتی‌های حاد تنفسی» سوزش‌های شدید در ناحیه نای» بیماری سیستم عصبی. سرطان ریه و پوست (رحمانیان و 
صفری. ۱۳۹۹) توسط سازمان حفاظت محیط‌زیست در گروه آلاینده‌های دارای تقدم قرار دارد (ززولی و همکاران؛ 
۱ ) این عنصر یکی از عناصر کمیاب است که از دو منبع طبیعی (هوازدگی سنگ‌ها و خاک) و انسان‌زاد (صنایع» 
کشاورزی و ploy‏ نقلیه) به محیط‌زیست وارد می‌شود (کمپل و QU PSU‏ 

مطالعات. غلظت بالای نیکل در خاک‌های تشکیل شده بر روی مواد مادری بازیک و فوق بازیک را نشان داده 
است (سالوادوربلانس" و همکاران, Y V‏ دآمیکو " و همکاران: ۲۰۰۸). همچنین مطالعه‌ی غلظت نیکل در خاک‌های 
مناطق شمال غرب Ol pl‏ نشان oala‏ است که بیش از AV‏ درصد از خاک‌های این مناطق دارای غلظت GUL‏ این عنصر 
نسبت به استانداردهای تعیین شده در ایران هستند و بیشترین مقدار JSS‏ در خاک‌های تشکیل شده روی مواد مادری 
فوق بازیک با متوسط VW‏ میلی‌گرم بر کیلو گرم مشاهده شد (یوسفی فر و همکاران, ۱۳۹۸). محققین نشان داده‌اند 
غلظت نیکل در خاک‌های حاصل از سنگ‌های فوق بازیک تا ۲۰۰۰ میلی‌گرم بر کیلوگرم نیز می‌رسد (الکساندں 
۶ جدول (۱) ميانگین غلظت نیکل را در برخی از انواع سنگ‌ها نشان می‌دهد. 


1 Islam et al. 

2 Kasassi et al. 

3 Caillaud et al. 

4 Sari et al. 

5 Cempel & Nikel 
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7 D'Amico 
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حدول ۱- میانگین غلظت نیکل (میلی گرم بر کیلوگرم) در سنگ‌های مختلف (آلووی. ۱۹۹۰) 


فلز 


سنگ‌های آذرین 
m‏ 


jo: 


گرانیت 


۰/۵ 


سنگ‌های رسوبی 
سنگ آهک ماسه‌سنگ 


q ۷ 


yyy 


آلاینده‌ها در جهت جریان وزش باد حرکت می‌کنند و در هوا معلق مانده و به‌مرورزمان بر روی SLE‏ و OLALS‏ 
نشست می‌کنند؛ بنابراین طبیعی است که میزان آلودگی pole‏ سنگین در جهت وزش باد به‌ویژه باد غالب بیشتر باشد 
(آگولینی | و همکاران ۲۰۱۳). 
احتمالاً مصرف آفت‌کش‌ها و کودهای شیمیایی به‌ویژه کودهای فسفاته در مناطق کشاورزی در بالا رفتن غلظت 
نیکل کل و قابل‌جذب در این اراضی مؤثر بوده است. هرساله بیش از ۳۰ میلیون تن کودهای فسفاته در سرتاسر 
جهان مورد استفاده فرار می‌گیرد که بیش از ۹٩‏ درصد آن‌ها از سنگ‌های فسفاته حاصل می‌شوند. مشکل اساسی که 
uan‏ گت wig‏ کرای قات در Os Sat tet cesa eS eli‏ این E tl‏ 
همچون نیکل است ( لامبرت و همکاران. ۲۰۰۷). جدول Y‏ غلظت نیکل را در چند نوع کود فسفاته OUS‏ می‌دهد. 


حدول ۲- میانگین غلظت نیکل در چند نوع کود فسفاته (لامبرت " و همکاران. ۲۰۰۷) 


نوع کود 
مونوآمونیوم فسفات 
دی آمو نیوم فسفات 
سوپر فسفات des‏ 


سنگ فسفات 


غلظت نیکل (میلی گرم بر کیلوگرم) 
۲۲/۲ 
evr‏ 
۱۷/۳ 


0/4 


از سوی دیگر استفاده از فاضلاب در آبیاری اراضی کشاورزی بخصوص فاضلاب‌های صنعتی نیز در بالا رفتن 
غلظت نیکل اراضی کشاورزی بسیار موثر است (آتش پر و رضا پور (WA‏ 
غلظت JSS‏ نیز مانند سایر متغیرهای محیطی تحت تأثیر عوامل متعددی است که به‌صورت مکانی و زمانی تغییر 


می‌کنند و این امر شناخت گسترش و توزیع این متغیرهای محیطی را با پیچیدگی زیادی روبه‌رو کرده است (وو "و 


1 Ugolini et al. 
2 Lambert et al. 
3 Wuet al. 
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همکاران. ۲۰۲۰). لذا انجام توصیه‌های مدیریتی و تدوین برنامه‌های مرتبط با این متغیرها در محیط‌زیست مستلزم 
تعیین الگوی تغییرات مکانی آن‌ها است. ازاین‌رو به دلیل مشکلات مربوط به نمونه‌برداری» استفاده از روش درون‌یابی 
بسیار مفید است (حلیمی و همکاران» ۱۳۹۲). در این روش‌ها می‌توان با جمع‌آوری داده‌های کمتر نقشه‌های 
دقیق‌تری را a‏ نمود که می‌تواند به مدیریت صحیح و اصولی منجر گردد. از میان همه‌ی روش‌های درون‌یابی» 
روش زمین‌آمار محبوبیت خاصی دارد. این روش یکی از دقیق‌ترین روش‌هایی است که علاوه بر توصیف تغییرات 
مکانی و Shy‏ داده‌هاء قادر به تهیه نقشه‌های کمی پراکنش متغیرهای محیطی با حداقل واریانس ممکن است 
(گوائرتز ؛ ۲۰۰۰؛ موکانا و کوئیکی T‏ ۲۰۰۸)؛ اما تاکنون از بین تمام روش‌های زمین‌آمار روش کریجینگ معمولی به 
دلیل cles‏ کم‌تر و تولید حداقل واریانس ممکن» محبوبیت بیشتری را نشان داده است و در تعیین توزیع مکانی 
متغیرهای محیطی پتانسیل بالایی دارد. در زمینه استفاده از روش‌های زمین‌آمار در بررسی توزیع مکانی pole‏ در 
خاک مطالعات متعددی در Us‏ (مانوز Leb‏ و همکاران» ۲۰۲۰؛ FOS‏ و همکاران»۲۰۱۹؛ بهونیا* و همکاران ۲۰۱۸؛ 
کائو" و همکاران» ۲۰۱۷) و ایران (صدر و موحدی راد ۱۶۰۰ صدر و موحدی caly‏ ۱۳۹۹ محمدی. WV‏ صدر و 
افیونی» Y‏ استواری و همکاران ۱۳۹۶؛ حداکرمی و همکاران, ۱۳۹۳؛ میرزایی و iab Saa‏ ۱۳۹۲؛ عبداللهی و 
همکاران» ۱۳۹۱؛ شهبازی و همکاران. ۱۳۹۰ ) انجام گرفته است. 

علی‌رغم اهمیت GUL‏ وجود اطلاعات کافی و کامل از وضعیت توزیع عناصر سنگین در استان‌های کشوره جهعت 
اتخاذ تصمیمات درست مدیریتی در پروژه‌های ملی و منطقه‌ای» تنها در تعداد معدودی از استان‌های کشور پراکنش 
فلزات سنگین به‌صورت نقشه‌های کاربردی ارائه شده است و در بیشتر استان‌هاء ازجمله استان قم اطلاعات بسیار 
محدودی در مورد مقدار کل فلزات سنگین در خاک وجود دارد. بنابراین» به دلیل تمرکز فعالیت‌های صنعتی و 
کشاورزی. تعیین وضعیت پراکنش pole‏ سنگین در خاک‌های سطحی این منطقه بسیار حائز اهمیت است. لذا 
پژوهش حاضر با اهداف ذیل انجام گرفت. ۱- تعیین توزیع مکانی نیکل در بخشی از اراضی استان قم. ۲- تعیین 
مکان‌های آلوده به نیکل در مقایسه با استانداردهای موجود در مورد آن. ۳- بررسی مهم‌ترین منابع ورود این فلز به 
مناطق آلوده. 


1 Goovaerts 

2 Moukana & Koike 
3 Mufioz-Nájera et al. 
4 Zhen et al. 

5 Bhunia et al. 

6 Cao et al 


سال يازدهم توزیع منطقه‌ای نیکل در SSE‏ سطحی .... YA‏ 


۲- مواد و روش‌ها 

۱-۲- موقعیت جغرافیایی منطقه مطالعاتی 

منطقه مطالعاتی بخش‌هایی از اراضی کشاورزی» شهری-صنعتی و بایر استان قم به مرکزیت شهر قم» به وسعت 
۶ کیلومترمربع بین طول‌های جغرافیایی 4۷ ۵۰° تا ۸ 019 شرقی و عرض‌های جغرافیایی 40 Yt?‏ تا ۳۰ ۳۹۶ 
شمالی قرار دارد. شکل LIL)‏ تصوير ماهواره‌ای منطقه مطالعاتی را نشان می‌دهد. میانگین ارتفاع منطقه از سطح دریا 
۰ متر و میانگین بارندگی در منطقه مطالعاتی بر اساس سالنامه آماری استان در یک دوره ۲۶ ساله در حدود ۱۳۸ 
میلی‌متر با حداکثر ۲۰ میلی‌متر در سال است. جهت وزش بادهای غالب در منطقه از سمت غرب است (شکل ۱- 
ب). به‌طورکلی تشکیلات زمین‌شناختی غالب منطقه از سازندهای جوان دوران چهارم و سنگ‌های ترشیاری و 
کواترنر تشکیل گردیده که عمدتاً شامل رسوبات آبرفتی و تپه‌های ماسه‌ای دوران چهارم و سنگ‌های رسوبی و 
آتشفشانی کمی دگرگون شده دوران سوم است (شکل ۲). این منطقه بر روی نقشه سازمان جغرافیایی کشور با 
مقیاس ۱:۵۰۰۰۰ بانام قم باشماره 11 1۱۵۹ و قمرود با شماره TIE‏ 1۲۵۹ مشخص شده است 


http://www.giscoders.com) 


n الف‎ 


شکل ۱- الف. تصویر ماهواره‌ای منطقه مطالعاتی حهت تهیه نقشه کاربری اراضی TM)‏ 4-7 -2) و ب. گلباد 
منطقه مطالعاتی 


YY 


۱۱ 
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پادگانه‌های آبرفتی کهن 

پادگانه‌های آبرفتی جوان 

ej 

گدازه متوسط بازیک شيشه PERS‏ آتشفشانی و اسیدی 
کفه رسی رس و سیلت 

کنگلومرا با درون لایه‌های ماسه سنگی و رسی 

آبر فت wa>‏ قلوه و ماسه 


تراورتن 
سنگ‌های اتشفشانی base‏ توف 
مارن قرمز همراه با لایه‌های ماسه سنگی» شیلی. کنگلومرایی و گچی 


شکل ۲- نقشه زمین‌شناسی منطقه به تفکیک محدوده منطقه مطالعاتی و با مقیاس ۲۰۰۰۰ع:۱ 


-Y-Y‏ نمونه‌برداری 


در این مطالعه به خاک سطحی (عمق ۱۰-۰ سانتی‌متر) توجه شده است. نمونه‌برداری در Y*‏ نقطه بر روی 


شبکه‌ای با فواصل حدود ۱/۵*۱/۵ کیلومتر در اراضی کشاورزی» شهری-صنعتی و حدود VY‏ کیلومتر در اراضی 
بایر انجام شد. موقعیت جغرافیایی نمونه‌ها توسط دستگاه GPS‏ تعیین و کاربری محل نمونه‌برداری ثبت شد. شکل ۳ 


موقعیت نقاط نمونه‌برداری شده را با تفکیک نوع کاربری نشان می‌دهد. در این مطالعه ۱۰۳ نمونه در اراضی با 


کاربری کشاورزی و ۳۱ و ۷۵ نمونه به ترتیب در اراضی با کاربری شهری-صنعتی و بایر قرار گرفت. 


0 25006 


شکل ۳- Cond ge‏ نقاط نمونه‌برداری شده در منطته مطالعاتی الف. با تفکیک نوع کاربری اراضی منطقه ب. با 


توجه به موقعیت جغرافیایی نقطه 


سال یازدهم توزیع منطقه‌ای نیکل در خاک‌های سطحی .... ۳۳ 


۳-۲- اندازه گیری غلظت JSS‏ کل و قابل جذب در نمونه‌ها 

یک گرم از نمونه‌های خاک هوا خشک توزین و در ارلن مایر ۲۵۰ سی‌سی. ۲۰ میلی‌لیتر اسید نیتریک ۵ نرمال به 
آن‌ها اضافه شد. محتویات ارلن به مدت ۱۲ ساعت به حال خود رها شدند و سپس به مدت ۱۰ دقیقه جوشانده 
شدند. بعد از سرد شدن نمونه‌ها توسط کاغذ صافی واتمن ۶۲ در بالن ژوژه‌های ۲۵ میلی‌لیتری صاف گردیده و با 
آب مقطر به حجم رسانده شد. درنهایت غلظت کل JSS‏ توسط دستگاه جذب اتمی مدل Perkin Elmer‏ قرائت 
گردید(آجای و کمسون" ۱۹۸۳) 

شکل قابل جذب JSS‏ به‌وسیله محلول DTPA‏ ۰/۰۰۵ نرمال حاوی کلرید کلسیم ۰/۰۱ نرمال و تری اتانول 
آمین ۰/۱ نرمال اندازه‌گیری شد. pH‏ محلول DTPA‏ با استفاده از HCI‏ نرمال در ۷/۳ تنظیم گردید. به این صورت که 
تعلیق با نسبت Y‏ به Y‏ خاک به محلول DTPA‏ تهیه و بعد از ۲ ساعت همزدن نمونه‌ها با سرعت ۱۲۰ دور در BBS‏ 
از کاغذ صافی واتمن EY‏ عبور داده شد. سپس غلظت شکل قابل جذب نیکل توسط دستگاه جذب اتمی اندازه‌گیری 
شد (لینزی و نارول  (VA‏ 

در این مطالعه برای صحت‌سنجی نتایج اندازه‌گیری غلظت نیکل در خاک. یک نمونه استاندارد از موسسه 
تحقیقات ملی و فناوری آمریکا (MST)‏ به نام 2709 # ag San Joaquin‏ گردید. غلظت نیکل نمونه مذکور t AA)‏ 
۵ گزارش شد و به همان طریق شکل کل و قابل جذب JSS‏ نیز اندازه‌گیری شد. بازیابی مقدار فلزات آندازه‌گیری 
شده از نمونه استاندارد برای als‏ حاک‌های موردمطالعه بیش از ٩۰‏ درصد بود. 

بر اساس استاندارد WHO‏ غلظت نیکل را در خاک‌های آلوده ٩‏ و آلودگی متوسط و کم را به ترتیب ٦۳ر‏ 
میلی‌گرم بر کیلوگرم اعلام شده است (مکدونالد " و همکاران:۲۰۰۰). سازمان محیط‌زیست ایران مقادیر ۱۱۰ و ۵۰ 
میلی گرم بر کیلوگرم را به‌عنوان مقادیر استاندارد نیکل برای کاربری‌های کشاورزی و حفاظت محیطزیست گزارش 
کرده است (سازمان حفاظت محیطزیست ایران). جدول(۳) حداکثر مقدار مجاز نیکل را در خاک بر اساس 


جدول ۳- حداکثر مقدار مجاز JSS‏ در خاک برحسب میلی‌گرم بر کیلوگرم(پایس و بنتون جونز“ 1۹۹۷) 
عنصر کمیته اقتصادی اروپا هلند آلمان شمال شرقی امریکا 


Yo 0۰ TE Yo نیکل کل‎ 


1 Ajayi & Kamson 

2 Lindsay & Norvell 

3 MacDonald et al. 

4 Pais & Benton Jones 


ve‏ جغرافیا و مخاطرات محبطی شمارة اول 


حداکثر غلظت قابل‌قبول JSS‏ مطابق با استاندارد کشورهای مختلف با توجه به سمیت آن برای OLALS‏ در 


خاک‌های کشاورزی در جدول É‏ نشان داده شده است. 


جدول - حداکثر غلظت قابل قبول JSS‏ برحسب میلی گرم بر کیلوگرم در کشورهای مختلف با توجه به سمیت 
آن‌ها برای گیاه در خاک‌های کشاورزی (پایس و بنتون جونز؛ ۱1۹۹۷) 
فلز استرالیا کانادا لهستان a5‏ انگلستان آلمان 


Yes 5۰ £e ۱۰۰ m m نیکل‎ 


۲-- تجزیه و توصیف داده‌ها 

به‌منظور توصیف مشاهدات و نتایج آزمایش‌ها و دست‌یابی به خلاصه‌ای از اطلاعات آماری توصیف داده‌ها در 
lao‏ نرم‌افزار 16 SPSS‏ انجام گرفت. برای بررسی تست نرمال بودن داده‌ها در سطح اطمینان ۵ درصد از آمار 
کلم وگروف-اسمیرنوف استفاده شد. جهت ترسیم و تعیین مدل تغییرنما از برنامه رایانه‌ای Variowin‏ استفاده گردید. 
در این مطالعه L‏ تست a Cross-Validation‏ کمک نرم‌افزارء ٥۳‏ بهترین مدل تغییرنما که کمترین خطای 
تخمین را ایجاد نماید. انتخاب و میان‌یابی با توجه به آن‌ها به روش کریجینگ معمولی با همان برنامه رایانه‌ای 
صورت گرفت. در این مطالعه صحت الگوی برازش esl‏ شده به تغییرنما توسط میانگین مربعات خطای تخمین 
(MSE)‏ و میانگین خطای تخمین (ME)‏ سنجیده شد. به‌این‌ترتیب که این دو خصوصیت برای هرکدام از 
تغییرنماهای برازش داده شده به دست dal‏ تغییر نمای مناسب بر اساس مقدار این دو خصوصیت انتخاب گردید. 
مقدار پارامتر ME‏ درصورتی که صفر یا نزدیک صفر باشد نشانگر نااریب بودن تخمین است (حلیمی وهمکاران, 
۲ و حداقل بودن MSE‏ نشانگر بالا بودن دقت تخمین است (صدر و موحدی راد ۱۳۹۹). 

نقشه کریجینگ هم توسط برنامه رایانه‌ای 16 Surfer‏ تهیه شد. در این مطالعه کریجینگ نقطه‌ای HS‏ در خاک با 
ابعاد بلوک‌های 7۱۰۰۰ ۱۰۰۰ متر و با استفاده از الگوی تغییرنمای همه جهت و در نظر گرفتن ناهمسانگردی انجام 
گرفت. دامنه همبستگی تغییرنما در جهات مختلف احیه جستجو به شکل بیضی با قطر بزرگ با زاویه ٩۰‏ درجه از 
محور عمودی در نظر گرفته شد. در داخل aol‏ جستجو حداکثر ۸ و حداقل daz ٤‏ به‌منظور میان‌یابی هر نقطه بکار 


رفته است. 


1 Pais & Benton Jones 
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۳- نتایج و بحث 

۱-۳- نتایج آماری 

میانگین غلظت JSS‏ کل و قابل‌جذب در منطقه مطالعاتی به ترتیب VV‏ و ۰/۵ میلی‌گرم بر کیلوگرم بوده و 
ضریب تغییرات JSS‏ کل ۷۳/۹ درصد با حداقل ۰/۲ و حداکثر ۵۰ میلی گرم بر کیلوگرم و ضریب تغییرات SSE‏ 
قابل جذب Ve‏ درصد با حداقل ۰/۰۳ و حداکثر ۱/۸ میلی گرم بر کیلوگرم است COS pte)‏ در پژوهشی در خاک‌های 
استان اصفهان غلظت نیکل کل به‌طور متوسط ۰ میلی‌گرم بر کیلوگرم گزارش شده است. 


جدول ۵- خلاصه‌ای از وضعیت آماری غلظت کل و pU‏ جذب ISS‏ در منطقه مطالعاتی به تفکیک کاربری 


n i‏ . ۳ ضریب 
متغیر mg kg"‏ کاربری تعداد oly‏ واریانس حداکثر حداقل چولگی . 
تغییرات C1)‏ 
کشاورزی ۱۳ VAA Yo Vo Os VN V/A‏ 
شهری ۳۵ VA ۱۸ ۳۹/۶ YA/A NY‏ 2 
غلظت نیکل کل 
ovo ۰/۹ Y Yo ۳۳/۳ ۸/۵ VY E‏ 
کل ۳۹ VY/A YN ۰/۲ Os 0/۱ MA!‏ 
کشاورزی vw it A VÉ ۸ A ev‏ 
شهری vw A ^ VE TA A Yo‏ 
غلظت نیکل 
VE hy VA Y n ۷ ub‏ ۱ 
قابل جذب 
کل ۳۹ 0/ IN LY Y^ ۰/١‏ 1 


ضریب تغییرات (CV)‏ نشانگر درجه تغییرپذیری یک ویژگی در خاک است. درصورتی که % 0۷520 باشد 
نشان‌دهنده تغییریذیری اندک» 50% CV«‏ 2196€ تغییریذیری متوسط و 100% CV«‏ 5096€ تغییریذیری بالا را 
OUS‏ می‌دهد (کریمی نژاد ' و OL as‏ ۲۰۱۵). در این پژوهش ضریب تغییرات نیکل در هر دو شکل موردمطالعه 
بالاتر از :۵ درصد تعیین شده است که بیانگر وجود تغییرات obj‏ این متغیر در خاک منطقه موردمطالعه است. ضریب 
تغییرات JSS‏ در اطراف کارخانه سیمان در پاسوج کمتر از ۵۰ درصد گزارش شده است که نشان‌دهنده تغیبرات کم 
غلظت این عنصر در خاک است؛ درحالی که غلظت کادمیوم در این مناطق ضریب تغییرات UL‏ ۱۰۰ درصد را نشان 
داده است که بیانگر تغییرات بسیار obj‏ این عنصر در خاک منطقه مطالعاتی است (رحمانیان ۱۳۹۹). گزارش شده 
است که تغییرات زیاد یک متغیر می‌تواند نقش عوامل مدیریتی نسبت به عوامل موروثی در منشأ آن متغیر را بالاتر 


ببرد (فروغی‌فر و همکاران. OTe‏ 


1 Karimi Nezhad et al. 


YA‏ جغرافیا و مخاطرات محیطی شمارة اول 


تفاوت بین حداکثر نیکل کل (۵۰ میلی گرم بر کیلوگرم) و حداکثر نیکل قابل‌جذب VN)‏ میلی‌گرم بر کیلوگرم) 
به‌خوبی نمایانگر این وضعیت است که بخش قابل‌توجهی از نیکل در خاک‌های منطقه به حالت تثبیت شده وجود 
داشته و از دسترس گیاه دور بوده است و بخش بسیار کمی از نیکل تثبیت شده در بخش‌های مختلف خاک توسط 
DTPA‏ زاد شده Ces‏ تیمورزاده و همکاران AYMA)‏ در بررسی مکانی نیکل در شهرستان اسفراین غلظت نیکل کل 
اعلام شده است (لیو" و همکاران» ۲۰۰۳). 

پر اساس حدود مجاز اعلام شده توسط کشورهای مختلف GVO)‏ ۰ میلی گرم بر کیلوگرم» خاک‌های منطقه 
کشاورزی توجه ویژه‌ای را می‌طلبد تا از UL‏ رفتن غلظت JSS‏ از حدود مجاز توصیه‌شده جلوگیری شود. نتایج این 
مطالعه نشان داد که کمی بیش از ۶/۵ درصد نقاط نمونه‌برداری شده داری نیکل کل بالاتر از ۲۵ هستند که بر اساس 
حدود مجاز کمیته اقتصادی ارویا و شمال شرقی آمریکا آلوده محسوب می‌شوند. 
آزمون کولموگرفسمیرنف هم این مطلب را تأیید کرد؛ اما ازآنجایی که برای انجام محاسبات زمین‌آماری فرض 
نرمال بودن توزیع داده مطرح است. لذا لازم بود که توزیع نرمال گردد. در میان ۹ داده موحود تعداد ٩‏ داده 
به‌صورت نامتعادل eal)‏ پرت) بودند و باعث افزایش چولگی و انحراف توزیع از حالت نرمال گردیدند و به دلیل 
اهمیت استفاده از داده‌های اصلی. ترجیح داده شد که در محاسبات بعدی این داده‌ها از مجموع کل داده‌ها Sie‏ 
گردیده و در محاسبات بعدی از ۲۰۰ داده استفاده گردد. SS‏ این «So‏ نیز ضروری است که داده‌های حذف‌شده ۶/۳ 
درصد کل داده‌ها را تشکیل دادند. وجود داده‌های پرت از طریق نمودار جعبهای (شکل *ج) و وجود شکاف و 
احتلاف aL;‏ در نمودار توزیع فراوانی(شکل ۶-الف) بهخوبی مشخص است. نمودار توزیع فراوانی نرمال شده پس از 
حذف داده‌های co Ste‏ در شکل ۶ب OUS‏ داده شده است. 
نیکل کل را دارا می‌باشند. کاربری شهری-صنعتی از نظر غلظت نیکل کلبه‌طور معنی‌داری کمتر از کاربری کشاورزی 


بوده اما اراضی pl‏ دارای اختلاف معنی‌دار با دوکاربری دیگر نداشت (شکل٤-ج)‏ 


1 Liu et al 
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شکل ۶- الف. توزیع فراوانی نیکل کل ب. توزیع فراوانی نیکل کل پس از حذف داده‌های پرت ج. نمودار 
جعبه‌ای نیکل کل 


با توجه به شکل ILO‏ وضعیت پراکنش نیکل قابل‌جذب در منطقه از تابع نرمال پیروی کرده و چجولگی ۰/۸ و 
آزمون کولموگرفسمیرنف هم این مطلب را تأیید می‌کند. شکل (۵-ب) نمودار جعبه‌ای و مقایسه میانگین غلظت 
نیکل bE‏ جذب در کاربری‌های مختلف را OLS‏ می‌دهد. اراضی کشاورزی و شهری-صنعتی با میانگین ۰/۱ 
میلی‌گرم بر کیلوگرم از لحاظ JSS‏ قابل جذب در یک سطح قرار دارند. با توجه به شکل (۵-ب) میانگین SSS‏ 
قابل جذب در اراضی کشاورزی و شهری-صنعتی دارای اختلاف معنی‌دار نبوده ولی این دو کاربری در سطح 


معنی‌داری میانگین نیکل قابلجذب بیشتری نسبت به اراضی ab‏ دارند. 


YA‏ جغرافیا و مخاطرات محیطی شمارة اول 
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نیکل قابل جذب (میلی‌گرم بر کیلو گرم) شهری yl‏ کشاورزی 
ل ud whe cda:‏ کل ie UU.‏ ت yd‏ چان ل Mi‏ جات 


۲-۳- محاسبه و الگوسازی تغییرنما 

به‌طورمعمول در عملیات زمینآماری با فرض عدم وجود روند. اقدام به محاسبه تغییرنما می‌شود. گاهی برخی از 
تغییرنماها با افزایش فاصله در محدوده موردمطالعه همواره افزايش می‌یابند. این تغییرنماها می‌توانند دلالت به وجود 
روند" در محدوده موردمطالعه داشته باشند (محمدی» ۱۳۸۵). لیکن در برخی از شرایط فرض عدم وجود روند 
معتبر نیست و حضور روند در مقادیر داده‌ها باعث به‌هم‌ریختگی ساختار تغییرنما تجربی می گردد (مدنی؛ ۱۳۷۶). 
عموماً شناسایی روند در آلودگی با ترسیم تغییرات فاصله در مقابل غلظت به دست می‌آید. در صورت وجود روند 
در یک بعد خاص, مقادیر پارامتر مورد نظر دارای توزیع یکنواختی نبوده بلکه از مدل خاصی تبعیت می‌کنند.در این 
تحقیق نمودار پراکنش JSS‏ در دو فرم کل و قابل‌جذب در مقابل طول و عرض جغرافیایی منطقه مطالعاتی به‌عنوان 
معیاری از فاصله ترسیم شد (اشکال ٦‏ و ۷) و روند مشخصی در توزیع غلظت عناصر مشاهده نگردید. 


غلظت نیکل کل 
che)‏ گرم بر کیلو گرم) 


45000 عرض جغرافیایی (متر) 3815000 51500 طول جغرافیایی(متر) 47000 
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3848000 عرض جغرافیابی O)‏ 3113000 


(8 طول جنرافياي مت‎ 3M 


شکل ۷- نمودار پراکنش غلظت نیکل قابل جذب در الف. عرض جغرافیایی و ب. طول جغرافیایی 


۳۹ 


در مطالعه حاضر با توجه به طول منطقه مطالعاتی که ۳۱ کیلومتر بود. حداکثر فاصله مجاز در رسم تغییرنما VV‏ 
کیلومتر انتخاب شد. فاصله از ۰/۵ کیلومتر تا ٤‏ کیلومتر که حداقل و حداکثر فاصله بین نمونه‌ها هستند به‌عنوان گام 
برای مطالعه تغییرات مکانی موردبررسی قرار گرفت. تحمل به فاصله نصف فاصله در نظر گرفته شد. در این مطالعه 
با استفاده از سعی و خطا مناسب‌ترین الگوهای تغییرنمای متغیرهای موردمطالعه تعیین گردید و سپس توسط آنالیز 
کنترل اعتبا اعتبار هر تغییرنما به دست آمد و از میان حداقل ۵ و حداکثر ۱۸ الگی یک الگو به‌عنوان بهترین الگو 
جهت کریجینگ انتخاب گردید. در این مطالعه جهت محاسبه معیارهای آماری از ۲۰۰ داده استفاده شد. این ۲۰۰ داده 


اصلی دارای مقدار مشخصی بودند. آنگاه در هر مرحله یک داده حذف و سپس مقدار آن توسط کریجینگ معمولی و 


با کمک الگوی تغییرنمای منتخب و سایر نمونه‌های اطراف OF‏ تخمین زده شد. برای مقایسه داده‌های واقعی با 
تخمین حاصل از کریجینگ» از ضریب همبستگی پیرسون استفاده شد که در سطح | درصد ھمبستگی نسبتاً بالایی را 
نشان داد و تأیید کننده دقت GUL‏ تخمین و اعتبار تغییرنمای منتخب است. جدول نتایج حاصل از مقایسه 


تخمین‌های حاصل را با داده‌های واقعی در نیکل کل و قابل‌جذب OUS‏ می‌دهد. 


حدول 7 مقایسه تخمین‌های حاصل از کریحینگ نیکل کل با داده‌های واقعی 


ویژگی 
نیکل کل 
یکل قابلجذب 
** - در سطح ۰/۰۱ معنی‌دار است. 


معیار 
ضریب پیرسون 


ضریب پیرسون 


ضریب همبستگی بین داده‌های واقعی و کربجینگ نقطه‌ای همان نقاط 
gem‏ 


(Joy 


مشخصات تغییرنمای منتخب مربوط به ویژگی‌های موردیررسی در جدول ۷ ارائه شده است. شدت و درجه 


وابستگی مکانی یک متغیر ناحیه‌ای را می‌توان از تقسیم واریانس اثر قطعه‌ای به واریانس کل (حد آستانه) به دست 


جغرافیا و مخاطرات محیطی 


نماز ة اول 


آورد. اگر نسبت مذکور کمتر از ۲۵ درصد باشد متغیر دارای کلاس وابستگی مکانی قوی و اگر نسبت بین ۲۵ L3‏ ۷۵ 
درصد ASL‏ آنگاه متغیر مورد نظر دارای کلاس وابستگی مکانی متوسط است. در مطالعه حاضر توزیع نیکل کل 
دارای وابستگی مکانی قوی و ISS‏ قابل جذب دارای وابستگی مکانی متوسط است (امینی " و همکاران, ۲۰۰۵). 
همان‌طور که در جدول V‏ مشاهده می‌شود مقدار MSE‏ و ME‏ نشان‌دهنده دقت مناسب تخمین است. همان‌طور که 


MSE‏ نشانگر بالا بودن دقت تخمین است. وسعت زیاد منطقه مطالعاتی و فواصل زیاد میان نمونه‌ها از عواملی است 


جدول ۷- پارامترهای تغییرنما متغیر نیکل کل در خاک 


" دامنه تأثیر | حد آستانه 
عنصر مدل ناهمسانگردی اثر قطعه‌ای (Co)‏ 
(Co C) | (Ao (m)‏ 
نیکل کل کروی :4 Y/AL‏ 0۸۰ ۲4/0 
نیکل قابل‌جذب | کروی .4 (JM 10.۰ T‏ 
+(|h‏ 
dh) E‏ 
à .‏ : 
۰ ۰ ۰ 
O ۳۷ » ۰‏ سس 
۳۴ 
۳۱ 
\A d‏ 
1۵ 
۱۳ 
af 7‏ 
۶ 


Ae Mer ۱۳۰۰۰ 
فاصله(متر)‎ 


شکل ۸- تغییرنمای الف. ني 


Yee — Pee £e. 


THY‏ میانیابی نیکل در خاک 


شکل ٩‏ و ۱۰ به ترتیب نقشه حاصل از کریجینگ نقطه‌ای ISS‏ کل و نقشه حاصل از کریجینگ نقطه‌ای JS‏ 
قابل جذب در خاک منطقه موردمطالعه می‌باشند. ضریب همبستگی بین JSS‏ کل و LSS‏ قابل جذب در این مطالعه 


درصد وابستگی مکانی 
Co + C/Co‏ 


\o 


۳۹ 


MSE ME 
Wít —N 
A TEE 
7۱0۱ 
۰/۱۲۶ | 
“ay | 
۰ 
۰ ۰ 
۶۲ ۰ 
۰ 
۰۰۳۱ | 


] Amini et al. 
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بسیار پایین )0/98( بود. گزارش شده که مقدار نیکل کل نمی‌تواند به‌عنوان معیار مناسبی برای تعیین میزان نیکل 
قابل جذب مطرح شود (بارکر و پیلبم » ۲۰۱۵). این موضوع را عدم انطباق نقاط داغ دو نقشه تأیید می‌کند. 


WwW 3 
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کا ر Go‏ به ترتیب حرم حضرت معصومه (س) و شهرک صنعتی محمودآباد است. 


شکل -٩‏ نقشه حاصل از کریجینگ نقطه‌ای نیکل کل (میلی گرم بر کیلوگرم) 


با توجه به نقشه کریجینگ نیکل کل و انطباق این نقشه با نقشه کاربری اراضی (شکل ۳-الف) بیشترین غلظت 
نیکل کل در قسمت‌های مرکزی» جنوب غربی و شمال غربی منطقه که کاربری کشاورزی» شهری صنعتی دارند» 
دیده می‌شود. توزیع ISS‏ کل از قسمت‌های شمال غربی با کاربری شهری صنعتی به سمت قسمت‌های مرکزی با 
کاربری کشاورزی UL‏ است و در قسمت‌های شرقی کمترین غلظت نیکل دیده می‌شود. نوع کاربری در افزايش 
غلظت نیکل کل تأثیر زیادی داشته است. 


1 Barker & Pilbeam 
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شکل ۱۰- نقشه حاصل از کریجینگ نقطه‌ای نیکل قابل جذب (میلی‌گرم بر GEIS‏ 


اگرچه غلظت نیکل در منطقه مطالعاتی بیش از حدود استاندارد جهانی نبود اما یکی از عوامل موثر در بالا رفتن 
غلظت نیکل در منطقه مطالعاتی را می‌توان به فعالیت‌های صنعتی و وجود کارخانه‌ها و کارگاه‌های کوچک و بزرگی 
که در سطح منطقه پراکنده هستند نسبت GAB ala‏ صنعت در گسترش و توسعه استان قم به‌مراتب بیش از بخحش 
کشاورزی بوده است. یکی از دلایل گسترش کارگاه‌های صنعتی در این استان» نزدیکی به شهر تهران و از طرف دیگر 
به دو استان pene‏ اصفهان و مرکزی است. در Ds‏ کسترگی بالای محدوده‌های صنعتی در منطفه بطالعاتی» شسهرک 
صنعتی محمودآباد به‌عنوان یکی از فعال‌ترین شهرک‌های صنعتی استان قم محسوب می‌شود که با وسعتی بالغ‌بر ۷۰۰ 
هکتار در بخش قمرود واقع شده است و ده‌ها کارگاه و کارخانه صنعتی را در 2 جای داده است (محل قرارگیری 
cil‏ کر eat‏ ای کف ا aed‏ ا هاو کر nici‏ در مدرد ما من تاه | 
می‌توان در این امر موثر دانست. بررسی‌های انجام شده روی SE‏ کارخانه سیمان یاسوج بر روی خاک‌های اراضی 
مرا E‏ کان دای کے ای کل ر سای ۱ مه ای ay‏ مس دراي EEE‏ 
استانداردهای جهانی USEPA‏ بود. به نظر محقق این سطح بالای نیکل نه‌تنها ناشی از انتشارات خود کارخانه بوده 
بلکه > ploy OS‏ نقلیه در اطراف کارخانه هم بی‌تأثیر نبوده است (رحمانیان» ۱۳۹۹). 

جهت وزش بادهای غربی نیز ازجمله عواملی است که درمجموع باعث افزايش نیکل حاصل از صنایع موجود 
as ail gs‏ شاک ad‏ اس در مطالعه po‏ دو et eas‏ از تیک کل وا و چاه ege‏ باه غاب فک 


yl cyl {=}‏ مشاهده می گردد. 


سال یازدهم توزیع منطقه‌ای JSS‏ در خاک‌های سطحی .... m‏ 


در این مطالعه نقاط داغ (غلظت بالای ۲۰ میلی گرم بر کیلوگرم غلظت نیکل) در اراضی با کاربری بایر قرارگرفته 
است و ازانجاکه سازندهای جنوب منطقه مطالعاتی با توجه به نقشه زمین‌شناسی (شکل ۲ از نوع آتشفشانی است 
به نظر می‌رسد این محدوده‌ی اندک که دارای غلظت JSS GUL‏ است» زبانه‌هایی است که تحت تأثیر نوع سازند 
مناطق جنوبی قرار دارد. احتمالاً مواد مادری یکی از عوامل اصلی افزایش ISS‏ به خاک به شمار می‌رود. زیرا حضور 
سنگ‌های آتشفشانی و همین‌طور شیل‌ها و سنگ‌های بازیک و فوق بازیک که غنی از نیکل هستند. در قسمت‌های 
شمالی و جنوبی منطقه می‌تواند عاملی در جهت افزایش JSS‏ به خاک‌های منطقه مطالعاتی محسوب شود. با توجه 
به نقشه کاربری اراضی (شکل ۲) بخشی از اراضی کشاورزی. در محدوده مناطق دارای نیکل قابل جذب بالاتر قرار 
می گیرند که تأثیر استفاده از نهاده‌های کشاورزی مانند کودها و سموم و آفت‌کش‌ها را در UL‏ رفتن غلظت این عنصر 
مورد تأیید قرار می‌دهد. برزین و همکاران VME)‏ در بررسی آلودگی برخی فلزات سنگین خاک‌های سطحی استان 
همدان گزارش کردند غلظت نیکل تحت JES‏ عوامل طبیعی مانند مواد مادری است. ولی فعالیت‌های کشاورزی به 
دلیل مصرف بیش ازحد کودهای دامی و شیمیایی می‌تواند باعث افزایش هر چه بیشتر این عناصر در خاک شود. 

البته لازم به ذکر است که غلظت نیکل قابل‌جذب بیش از اینکه در کل منطقه سمیت را نشان دهد نگرانی ناشی از 
کمبود را مطرح می‌کند. گزارش‌هایی وجود دارد که مشکلات شدید کمبود نیکل و کمبودهای پنهان این عنصر را در 
محصولات زراعی OUS‏ می‌دهد (آلووی (Y AU‏ شکل قابل جذب نیکل برای گیاه و تغذیه مناسب ol‏ به‌عنوان یک 
عنصر کم‌مصرف ضروری (براوون" و همکاران؛ ۱۹۸۷)» مطرح است. مشخص است که نقش اصلی نیکل در ترکیب 
آنزیم اوره‌از است و درنتیجه این عنصر برای کلیه گیاهان و به‌ویژه بقولات جهت انتقال نیتروژن از ریشه‌ها و انتقال 
به اندام هوایی و تولید ترکیبات اورئیدی لازم است (کاباتا پندیاس و پندیاس T‏ ۲۰۰۱) و کمبود این عنصر می‌تواند با 
افزایش فعالیت نیترات ردوکتانن غلظت نیترات گیاه را کاهش دهد (علی‌بخشی و خوشگفتارمنش* (Yo‏ 


é‏ جمع‌بندی 
به‌طورکلی نتایج نشان داد که غلظت JSS‏ در بخش‌های وسیعی از استان قم در محدوده سمیت قرار ندارد بلکه 
تگرانی ناشی از کمبود این عنصر را نیز ایجاد می‌کند؛ چراکه نیکل یک عنصر ضروری برای رشد LS‏ محسوب 
می‌شود. غلظت نیکل در غالب منطقه مطالعاتی در زير حد آلودگی و سمیت است؛ اما حدود ۵ درصد این اراضی 
دارای غلظت IS‏ بیش‌ازحد مجاز برخی کشورها مثل اتحادیه اروپا بوده و عدم توجه به کنترل ورود این عنصر به 
خاک می‌تواند مشکلاتی را برای ساکنین این مناطق ایجاد کند. نتایج نشان می‌دهد که تجمع صنایع در استان قم 


1 Alloway 

2 Brown et al. 

3 Kabata-Pendias & Pendias 

4 Alibakhshi & Khoshgoftarmanesh 
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توانسته غلظت نیکل را در SE‏ سطحی این استان تحت تأثیر قرار دهد چراکه غلظت این عنصر در کاربری‌های 
مختلف دارای تفاوت معنی‌دار بود. احتمالاً مواد مادری و سنگ‌های بازیک و فوق بازیک که غنی از نیکل هستند و 
استفاده از کودهای فسفاته و لجن فاضلاب باعث افزایش ISS‏ و جهت وزش باده ای غربی باعث انتقال نیکل در 
منطقه گردیده است. 


۵- تشکر و قدردانی 


با سپاس از دانشگاه ly‏ نور استان کرمان» مرکز رفسنجان. 
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